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We introduce a measurement-based method for verifying quantum discord of any bipartite quan- 
tum system. We show that if the reduced states of one subsystem, conditioned to all outcomes of 
an informationally complete POVM (IC-POVM) performed on the other subsystem, commute with 
one another, there is no quantum discord from the first subsystem to the second. In particular, this 
method enables us to verify quantum discord for continuous-variable systems, as IC-POVMs are 
readily available from homodyne and heterodyne measurements. Using this method, we show that 
the only Gaussian states with zero discord are product states; hence Gaussian states with Gaussian 
discord have nonzero quantum discord. 

PACS numbers: 42.50.Dv, 03.65.Ta, 03.65.Ud 



I. INTRODUCTION 

Quantum correlations play a central role as a re- 
source in quantum information processing and quantum 
communication tasks. Traditionally, entanglement was 
thought of as the unique form of quantum correlation 
and the reason why quantum computers can outperform 
classical computers. Yet there are tasks that are believed 
to be exponentially hard classically, which can be done 
efficiently using quantum computational models with lit- 
tle or no entanglement (T), ■ Quantum discord was in- 
troduced as a more general measure of quantum correla- 
tion for bipartite systems, with no classical analogue 
and it was suggested as a resource for certain quantum 
computation models [|[ . Discord has been generalized to 
continuous-variable systems to study quantum correla- 
tions in Gaussian states 0, @ and certain non-Gaussian 
states 0- 

Recently, some schemes have been proposed for wit- 
nessing nonvanishing quantum discord for discrete vari- 
able bipartite quantum states [H-[l3[ . Of particular prac- 
tical interest, however, is a general method for detecting 
nonvanishing discord in the joint state of any bipartite 
system. 

In this paper we introduce a measurement-based 
method for verifying quantum discord of any bipartite 
quantum state, without requiring any prior knowledge of 
the joint state. This method can be applied to both dis- 
crete and continuous-variable systems. We show that if 
the reduced states of one subsystem, B. conditioned to 
all the outcomes of an informationally complete POVM 
(IC-POVM) performed on the other subsystem, A, com- 
mute with one another, then the quantum discord from 
B to A is zero. Conversely, if they do not commute, the 
quantum discord from B to A is necessarily nonzero. A 
POVM is informationally complete if its outcome proba- 
bilities arc sufficient to determine uniquely the quantum 
state, i.e., to perform quantum state tomography flilfpoT ]. 



In other words, a bipartite state has zero discord from B 
to A if tomography on A leaves the eigenstates of the 
density operator of B unchanged. 

Quantum discord is defined as the difference between 
two classically equivalent measures for mutual informa- 
tion 0]. According to Bayes's rule for classically corre- 
lated probability distributions, the quantities I(A : B) = 
H(A) + H(B) - H(A,B), J{A\B) = H(A) - H(A\B), 
and J{B\A) = H(B) - H(B\A), where H denotes the 
Shannon entropy and H(A\B) = H{A, B) - H{B) is the 
conditional entropy, arc all equal; they are called the clas- 
sical mutual information. For a bipartite quantum sys- 
tem, the quantum mutual information is defined, in anal- 
ogy to I(A : B), by I(p AB ) = S(p A ) + S(p B ) - S{p AB ), 
where S{p) = — Tr[plog(p)] is the von Neumann entropy. 
A measurement-based, quantum version of the condi- 
tional entropy is S{ Uj }(A\B) = £V p j S(p A \ j ), where 
pj = Tr^sIIj], p A \j = Tr B [p AB Ilj]/pj, and the set 
{IP,-}, with ^2j^lj = I, makes up a POVM measure- 
ment on subsystem B. This conditional entropy depends 
on the choice of measurement; hence, the quantum ana- 
logue of J(A\B) is defined by minimizing over all possible 
measurements: J^{p AB ) = S(p A ) - M^yS^ji^lB). 

The quantum discord from B to A is then defined as 
the difference between these two ways of defining a quan- 
tum mutual information: 

D^(p AB ) = I(pab) - J^(pab) 

= S(p B ) - S(p AB ) + htf {nj} S {ni }(A\B) . 

(1) 

The quantum discord is zero if and only if the quantum 
state can be expressed in the form 



PAB = ^PjPj ® \j) (j 



(2) 



where {\j)} arc orthogonal states and < pj < 1 pjlllfjjl. 
For a quantum state with this form, local measurements 
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on B in the basis {\j ) (j |} leave the system unperturbed, 
and all the state information can be extracted without 
joint measurements. Notice that the state (0 is diagonal 
in a conditional product basis pointing from B to A, i.e., 
an orthogonal basis of the form {\fjk} <S> | ) } , where the 
states \fjk) are the eigenstates of pj. Hence, for a given 
quantum state, quantum discord can be directly verified 
by diagonalizing the joint density operator [l6j . 

Although there is no general method for minimizing 
over all possible measurements in order to calculate the 
conditional entropy of a state, this can be done by re- 
stricting to certain subspace of states and POVMs. The 
Gaussian quantum discord is defined as the quantum 
discord for two-mode Gaussian states where the evalua- 
tion of the conditional entropy is restricted to generalized 
Gaussian measurements [H, |g| . 

This paper is structured as follows. In the next section, 
we introduce our method for verifying quantum discord. 
In Sec. IIII1 we discuss the application of the method to 
continuous-variable systems. In Sec. IIV1 we show that 
for Gaussian states, only product states have zero dis- 
cord; hence states with nonzero Gaussian discord have 
nonvanishing quantum discord. 

II. METHOD 

Our method for verifying quantum discord is based on 
testing whether a quantum state can be expressed in the 
form @ of states with zero discord. This method is based 
on the following theorem. 

Theorem. For a bipartite system p AB , the necessary 
and sufficient condition for having zero discord from B 
to A, D^(pab) = 0, is that the states of subsystem B, 
PB\k = Tr A [M k p AB ]/ Tr AB [MkPAB] = Tr A [M k p AB }/pk, 
conditioned to the outcomes k of an IC-POVM on A 
(POVM elements {M k }), commute with one another, 
i.e., 

[pB\k, PB\k>] = , for any k and k' . (3) 
Proof. For a state having the zero-discord form @, 
E 3 P^A[M kPj ]\j)(j\ 

PB\k — ™ FT? 1 ' W 

immediately demonstrating that the states p B \ k are all 
diagonal in the basis and thus commute. 

For the converse, we assume the condition ([3]). For 
such a set of commuting conditional states {p B \ k }, there 
exists an orthonormal basis, that diagonalizes all 

the conditional states, p B \ k = ^kj \ j ) (j \- To say that 
the POVM elements {AI k } make up an IC-POVM is to 
say that they span the space of operators and thus there 
exist operators {N k } such that p = '^ k N k TT:[M kP \ for 
any density operator p. Applying this identity to the 
joint state gives 

Pab = N k Tr A [M kPAB ] = J2 ^kjPkN k ®\j) (i| . (5) 

k k,j 



Hence, p AB has the form ([2]), with pj = Xkjp k N k . 

□ 



Physically, what the proof says is that for a state of 
zero discord from B to A, the measurement of the POVM 
{M k <g> |j) (j |} extracts all information about the state 
p AB . From the perspective of the original definition of 
discord Q, one imagines extracting this information by 
first measuring B in the basis {|j)} and then measuring 
an IC-POVM on A. Our criterion for zero discord works 
from the opposite perspective by reversing the order of 
the measurements on A and B. 

In order to test whether an unknown quantum state 
has nonzero discord experimentally, based on this the- 
orem, one needs to measure an IC-POVM on subsys- 
tem A and determine, by state tomography for each 
outcome, the corresponding states of the subsystem B. 
This procedure continues until one of the commutation 
relations between conditional states of subsystem B is 
nonzero. For a quantum state in c?-dimensional Hilbert 
space, an IC-POVM has at least d 2 POVM elements. 
If some prior knowledge about the state in question is 
available, as is often the case in practice, quantum dis- 
cord can be checked with only a few IC-POVMs. Con- 
sider, for example, the maximally entangled state, 1^) = 
\j ) \j ) /Vd. Any two rank-one outcomes, |n) and 
|?7), on one of the subsystems, provided < | (n\rj) \ < 1 
(these could be outcomes from two distinct, nonorthogo- 
nal projective measurements), yield conditional states of 
the other subsystem do not commute. Also, as we show 
below, for Gaussian states only two different heterodyne 
outcomes arc sufficient to verify quantum discord. 



III. CONTINUOUS- VARIABLE SYSTEMS 

An interesting feature of this method is that it can 
be easily applied to continuous- variable systems, as com- 
plete sets of IC-POVMs are readily available from het- 
erodyne or homodyne detection. Two sets of measure- 
ments are required, one on each of the subsystems. In 
general, one needs to do state tomography to construct 
the quasiprobability distributions of subsystem B for all 
the states conditioned to outcomes of the measurement 
performed on subsystem A. Then the commutation rela- 
tions between states p B \ k , which are represented in terms 
of quasiprobabilities, must be calculated in order to verify 
discord. This can be done by using an appropriate rela- 
tion for calculating the commutation relations in terms 
of the reconstructed quasiprobabilities. For instance, if 
the Wigner functions of conditioned states of subsystem 
B, W B \ k (a) are available, the commutation relations be- 
tween corresponding density operators can be calculated 
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by using the Moyal Bracket [17| : 

W kk ,(a) = ^ f d 2 l3d 2 p'W Blk (a+±l3)W Blk ,(a+±l3>) 



(3 p'* -P'P* 



(0) 



Here W kk > (a) is the Wigner function for the operator 
~i[PB\ki PB\k']- If the states commute with each other 
then W kk ' (a) = for all a. Alternatively, the commuta- 
tion relations can be calculated using characteristic func- 
tions, 



x sin t 



(7) 



or in terms of any other quasiprobability distribu- 
tions [TJ. 

For states with zero discord, the eigenstates of the 
conditional density operator of B do not change while 
A is being fully determined from measurements of an 
IC-POVM. For continuous-variable systems, defined on 
an infinite-dimensional Hilbert space, the IC-POVM will 
have an infinite number of outcomes. In practice, only a 
finite number of measurement outcomes can be explored. 
For instance, in homodyne detection only a finite num- 
ber of phases are considered, and the phase space is sub- 
divided into a finite number of bins. This introduces 
errors in the state estimation and consequent uncertain- 
ties in verifying quantum discord. By having some prior 
knowledge about the state, however, the error can be es- 
timated, and it can be made arbitrarily small using a 
sufficiently large number of measurements. 



IV. GAUSSIAN STATES 

A special class of continuous- variable states consists of 
the Gaussian states, i.e., those states whose Wigner func- 
tion is a Gaussian function. Such states are uniquely 
characterized by the means and covariance matrix of 
their quadrature components, x and p. For two sys- 
tems, with modal annihilation operators a — x\ + ip\ 
and b = X2 + ip 2 , we define quadrature vectors for each 
system, xi = (x±,pi) and X2 = (x 2 ,p 2 ), and we define an 
overall quadrature vector x = (xi,x 2 ) = (xi,pi, x 2 ,p 2 ). 

The means of the quadrature components can be set 
to zero by locally displacing the two systems. Then the 
state is specified by its covariance matrix [Jj| 



(x T x) 



A C 
C T B 



(8) 



Using local unitary operations that preserve the Gaus- 
sian form of the states, the covariance matrix of a bipar- 
tite Gaussian state can be brought to a standard form in 



which A = diag(a, a), B = diag(fe, b), and C = diag(c, d), 
where a > and b > 0. This can be accomplishing 
by first applying local unitary rotations that diagonal- 
ize A and B, then using local squeezing operations to 
transform these diagonal blocks to A = diag(a, a) and 
B = diag(6, 6), and finally applying further local unitary 
rotations to diagonalize C. Notice that positivity of the 
density operator imposes the uncertainty-principle con- 
straint 12 Ol 



-O > 



n 











G 




J = 


5) 



(9) 



For a covariance matrix in standard form, this implies 
that a 2 > 1/16, b 2 > 1/16, ab > c 2 , and ab> d 2 , plus 
cubic and quartic constraints on a, 6, c, and d. 

It has been shown that Gaussian discord for a two- 
mode Gaussian state is zero if and only if C = @, 
6]. Here we show this condition is also necessary and 
sufficient for having zero discord. 

A zero-mean Gaussian state with the standard form of 
the covariance matrix has characteristic function x(k) = 
(e lkxT ) = e ~ XCTxT / 2 anc j Wigner function 

W(x lj x 2 ) = W r (x) 
1 

exp 



47r 2 V det i 



1 



^ 2 ^f(ab-c 2 ){ab-d 2 ) 



bx 2 + axn — 2cx\X2 bp 2 + api — 2dp\p2 \ 
^7-7 — ^7-7 — J 



2{ab-c 2 ) 



2{ab~d 2 ) 



(10) 



Suppose Alice makes a heterodyne measurement on 
subsystem A; i.e., she uses the IC-POVM whose POVM 
elements are the coherent states \(3)((3\. Let (3 = x[ +ip[ 
specify the outcomes of here measurement. Then the 
state Pb\xl[, conditioned on these outcomes, has Wigner 
function 



W B |xi(x 2 ) 



1 

N 



d 2 x 1 H^(x 1 ,x 2 )W xi (x 1 ) , (11) 



where W x ' (xi) = 2exp[— 2(xi — x[) 2 — 2(pi — p[) 2 ]/tt is 
the Wigner function of coherent state \/3) — \x[ + ip[) 
and N is a normalization factor. Integration yields 

W / i3|x' 1 (x 2 ) =-^ 7 exp^-ia; 2 /(a,6,c) + x 2 x' 1 g(a, b, c) 
X exp ( _ ^2/( a ) & ^) +P2p[g(a,b,d) 



(12) 



where 



f(a,b,z) = 
g(a,b,z) 



1 



ab — z 2 
4z 



b + A(ab- z 2 ) 



b + 4(ab-z 2 ) 
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and N' is a normalization factor. 

The peak of system £?'s conditional Wigncr func- 
tion (fT2|) . located at 

7 = 77 r ^ i +» 77 , tt Pi » (13) 
f{a,b,c) f{a,b,d) 

depends on the outcomes of the measurement on A, x[ 
and p'j, unless c = and d — 0. Consequently, the eigen- 
vectors of the conditional state Pb\x[- which are gener- 
ally displaced, squeezed number states, {D(^)S(Q \n)}, 
change based on the outcomes of the heterodyne mea- 
surement performed on subsystem A. Therefore, the bi- 
partite Gaussian state has nonvanishing discord unless 
c = and d = 0. Transforming back from the standard 
form to the general convariance matrix (j8J), one can say 
that a bipartite Gaussian state has zero discord if and 
only if C = 0, i.e., if and only if the state is a product 
state. This also implies that states with Gaussian discord 
(C 0) have nonzero quantum discord. 

These results show that quantum discord of Gaussian 
states can be verified using only two different heterodyne 
outcomes on one subsystem and finding (by tomography) 
the points in the phase space at which the corresponding 
conditional Wigncr functions attain their maximum val- 
ues. If those points do not coincide, the quantum discord 
is nonzero. This argument is also true for non-Gaussian 
states, as it is sufficient evidence that the corresponding 
operators do not commute. Note that discord exists even 
if only one of c and d is nonzero, so to uncover discord 
of Gaussian states with only two heterodyne outcomes, 
one should choose the outcomes to be different for both 
quadratures. 

V. CONCLUSION 

We have introduced a method for verifying quantum 
discord of any bipartite quantum system. The method 



is based on the fact that all the information in states 
with zero discord from subsystem B to subsystem A can 
be fully extracted by measurements that arc diagonal 
in a single basis of B. In order to verify discord, one 
needs to perform an IC-POVM on subsystem A and check 
whether the conditional density operators of the subsys- 
tem B commute, i.e. if they share the same eigenstates. 

It is worth mentioning that, in practice, one would cal- 
culate the commutation relations between the conditional 
states as an IC-POVM is being performed on A. Once 
one of these commutators is found to be nonzero, that 
confirms nonzero discord. This method can be simply 
applied on continuous-variable systems by using homo- 
dyne or heterodyne detection and calculating the com- 
mutation relations in terms of quasiprobability distribu- 
tions. We have shown that a bipartite Gaussian state has 
nonzero quantum discord if and only if it is not a prod- 
uct state, which is the same as the condition for having 
nonzero Gaussian discord. Moreover, we show that just 
two heterodyne outcomes are sufficient to verify quantum 
discord of Gaussian states. 
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